A little inflation in the early universe at the QCD phase transition 
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We explore a scenario that ahows for a strong first order phase-transition of QCD at non-negligible 
baryon number in the early universe and its possible cosmological observable consequences. The 
main assumption is a quasi-stable QCD-vacuum state that leads to a short period of inflation, con- 
sequently diluting the net baryon to photon ratio to it's today observed value. A strong mechanism 
for baryogenesis is needed to start out with a baryon asymmetry of order unity, e.g. as provided 
by AfHeck-Dine baryogenesis. The cosmological implications are direct effects on primordial density 
fluctuations up to dark matter mass scales of Mmax 1 — IOMq, change in the spectral slope up 
to mass scales of Mmax ~ (10* — 10** )Mq, production of primordial magnetic fields with initial 
strength up to Bqcd ~ 10^^ G and a gravitational wave spectrum with present day peak strain 
amplitude of up to hdi^peak) ~ 5 • 10^^^ around Vpeak ~ 4 • 10~**Hz. 
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The theory of quantum chromodynamics (QCD) 
predicts a phase transition from a quark-gluon plasma 
to a hadron gas in the early universe at a eritieal 
temperature Tqcd ~ 150-200 MeV [l|,ll. Only at low 
net baryon density lattice gauge theory indicates a rapid 
crossover from the quark-gluon-plasma to the hadronic 
phase. In the standard hot big bang scenario the baryon 
asymmetry is rjB ^ 10~^ — 10~^° already before the 
QCD phase transition and therefore the idea of a first 
order QCD phase transition in the early universe has 
been more or less abandoned. However, most of the 
QCD phase diagram is actually not well known. There 
has been recent progress in the attempt to include a 
finite baryon density on the lattice 0, |j] but in general 
one still has to rely on effective models (Hj to tackle 
the QCD phase diagram. However, a true first order 
phase transition is expected at finite baryon densities, 
as indicated by chiral effective models of QCD Q due to 
the melting of quark and/or gluon condensates and the 
phenomenon of color superconductivity Q. Therefore, 
wc would like to reopen the issue of a first order cos- 
mological phase transition by addressing whether there 
is a simple scenario in which the QCD phase transition 
at finite baryon densities can have consequences on 
cosmological scales. 

In this letter we demonstrate that the scenario of a 
little inflation at the QCD phase transition at high 
baryon densities is possible and not in contradiction to 
present cosmological observations. It has interesting 
cosmological implications though as it can directly affect 
primordial density fluctuations on dark matter mass 
scales below Mmax ~ 1 — IOMq, change the spectral 
slope up to mass scales of Mmax ^ (10^ — 10^)Mq due 
to the change of the global equation of state, produce 
primordial magnetic flelds that may be strong enough 
to seed the presently observed galactic and extragalactic 
magnetic fields and produce a spectrum of gravitational 
waves around a peak frequency of 4 • 10~^ Hz that rnay 
be observable via pulsar timing in the future [1, M,- 



Dark matter properties are also strongly affected as the 
annihilation cross section for cold dark matter has to be 
up to nine orders of magnitude lower to give the right 
amount of dark matter today, which can be probed at 
the LHC by detecting the neutralino with an unexpected 
low annihilation cross section, and thermal warm dark 
matter masses can be of the order of MeV without 
exceeding the decoupling degrees of freedom of the 
standard model. Such a cosmological phase transition 
would then bear more resemblance to the situation in 
heavy ion collisions or even the centre of neutron stars 
than to the standard QCD phase transition in the hot 
big bang scenario. Hence, the upcoming FAIR facility 
would actually be a probe for the physics of the early 
universe in this scenario. 

For a first order QCD phase transition in the early 
universe to be possible a nonvanishing baryochemical po- 
tential /is is necessary where hb/T ^ 0{1). The present 
day baryon asymmetry rjs = [tib ^ ^b)/''''^ h^-s been 
experimentally found to be 5.9 • 10~^° < Vb < 6.4 • 10~^° 
at 98% confidence by combining big bang nucleosyn- 
thesis, cosmic microwave background and large scale 
structure results [lo| . The number of baryons in a 
comoving volume is constant and can be estimated to 
be Nb ~ a^^BiTf ~ aj^BfT^ where the index i refers 
to the initial values when the vacuum energy starts 
to dominate over the radiation energy and / to the 
final values after reheating. Therefore the initial ratio 
of the chemical potential to the temperature can be 

higher by ^ ~ with = afja,. If the 

timescale for the decay of the false vacuum is short 
compared to the Hubble time then Ti ~ Tf and already 
for Q ~ 10'^ (corresponding to a little inflationary period 
with N 7 e-foldings) the initial baryon asymmetry r]Bi 
and iJ-i/Ti will be of order unity. Hence, the evolution 
of the early universe could pass then through the first 
order chiral phase transition of QCD. A well established 
mechanism for generating a high baryon number in the 
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early universe is the AfHcck-Dine baryogenesis The 
AfHcck-Dine mechanism produces in most cases far too 
much baryon number, thus either additional fields or 
more sophisticated coupling terms have to be introduced 
to reduce the initial baryon number production or it 
has to be reduced afterwards. For the latter case an 
obvious possibility would be a large entropy release 
that dilutes the baryon to photon ratio for example by 
an inflationary period (se e.g. ref. il3l)- Affleck-Dine 
baryogenesis can in fact produce 775 ~ where this 

is probably an upper bound [12 1. 

Wc note that the scenario proposed here has some 
similarities to thermal inflation as discussed by Lyth 
and Stewart [l3| as both are late time inflation periods 
in addition to ordinary inflation with a length of only 
about 10 e-foldings and may thus help to to resolve 
partly the moduli problem. In ref. [l^ the production of 
quark stars with masses of 10~^ — IOMq was proposed 
within a scenario similar to the one discussed here but 
at small baryon densities and without addressing the 
key consequences of such a second late time inflationary 
period. In [l5| it was recently proposed that a large 
lepton asymmetry could also result in a first order QCD 
phase transition in the early universe. 
An important issue of this approach is the stability of 
the barrier between the false and the true vacuum in the 
effective potential up to very low temperatures. This 
is indeed the case in chiral models of QCD including 
gluonic degrees of freedom in the form of a dilaton 
field in which case the barrier only vanishes in the 
T — >■ limit Csernai and Kapusta found in 

ref. 17| only small supercooling of about 1% below 



the critical temperature using values of the surface 
tension of about a ^ 50McV/fm^. The nuclcation rate 
r depends exponentially on the surface tension as well 
as on the free energy difference between both phases 
and it's ratio to the Hubble parameter T/H exhibits 
a maximum around ~ Tc/2. If F/iJ does not exceed 
one at this point the phase transition fails and we find 
that keeping the other parameters used in ref. [17[ the 
surface tension must exceed 450 MeV/im^ ^ 3.7T^. 
However, the precise value of the surface tension at the 
QCD phase transition at high densities is not known 
and has been a matter of debate, sec e.g. the discussion 
in ref. [l^ giving a possible range of cr = 50 to 150 
McV/fm^ without excluding even smaller or larger 
values. At very high densities calculations of the surface 
tension in the first order phase transition between color 
superconducting phases and nuclear matter arrive at 
surface tensions of up to 300 MeV/fm^ [l^.The value 
of the bag constant used by Csernai and Kapusta is at 
the upper end of values considered in the literature (i.e. 
B = (235MeV)^) and a reduction to the value found in 
the original paper of the MIT group by fits to hadron 
masses {B ~ (145McV)'*, see ref. [20|) also reduces the 
surface tension needed for nuclcation to fail to a value 



of 124 MeV/fm^. This of course only covers the initial 
failure to nucleate, but in general B and a will both be 
temperature dependent. After some supercooling (e.g. 7 
e-foldings at most) T/H must exceed one for inflation to 
end and the phase transition to occur. This could either 
take place due to a strong drop in the surface tension or 
even due to a complete vanishing of the barrier between 
the two phases in the effective potential. In the latter 
case the surface tension goes to zero and a spinodal 
decomposition takes place. This has been studied e.g. 
in [2l| for a bag like model. Strong sensitivities of 



nuclcation rates on the surface tension have been also 
found for high-density matter as encountered in the 
interior of neutron stars or in core-collapse supcrnovae 



22| so that nuclcation timescales can easily be in the 



range of /iS to the age of the universe. 
The equation of state has to fullfiU the usual condition 
€ + 3p < to enter an inflationary phase. In the Bag 
model this would be the case below a temperature 
Tinf = (30S/(p7r^))"^^^. In the linear-a-model or the 
NJF-model this occurs when the thermal contributions 
to the energy density become smaller than the vacuum 
contributions like the quark condensate (mqqq) « f'^m'i 
and thegluon condensate PQco/i^g) (G^^.G^'") « 4B. 
In ref. [23[ the idea of a " quench" in context of heavy ion 
collisions is discussed, i.e. the chiral field is trapped in 
a metastable minimum and supercools until the barrier 
in the effective potential disappears at zero temperature 
and the field "rolls down" to the true minimum. All 
in all a delayed chiral phase transition at high baryon 
densities can not be excluded for the early universe with 
our present poor knowledge of QCD at non-zero baryon 
densities. 

The majority of dark matter candidates is already 
chemically decoupled from the radiation fiuid at the 
QCD phase transition and thus do not participate in 
the reheating at the end of the inflationary period. 
Therefore the dark matter number density is diluted by 
the same factor 6^ as the net baryon number. Normally 
the dark matter mass enclosed inside the Hubble horizon 
is of the order of IO'^Mq at Tqcd ~ 170 MeV, so any 
influence on perturbations inside dark matter would not 
have any consequences on larger scales. An inflationary 
period at the QCD-phase transition can change this 
drastically, since the amount of dark matter enclosed 
inside the horizon must be larger by a factor 6^ initially 
to give the right amount of dark matter today. For a 
short inflationary period, as discussed here, there is an 
additional effect on perturbations that have physical 
wavenumbers kph < at the beginning of inflation. For 
general relativistic ideal fluid density perturbations the 
system of dynamical equations is closed by Einstein's 
i?o-cquation that reads (fc^^ -|- H)a = AttG{Sp + 3dp) in 
uniform expansion gauge (see e.g. [23 ]) where Sp and 5p 
are the sum of the density and pressure perturbations. 
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respectively and a is the perturbation of the lapse. The 
time derivative of the Hubble parameter is given via 
the second Friedmann equation 11 = — 47rG(e + p) = 

-47rG(|effi (^)^ + eur « (f )^ where the 

subscripts refer to matter and radiation with g = 3 to 
4, respectively, and the index % to the onset of inflation. 
Comparing this to the first Friedmann equation one finds 

that = 8^(ev- + effi(f)' + eM. (f)') meaning 

that the two scales differ by \H/H'^\'^/'^ ~ {^f''^ which 
would be irrelevant for a long infiationary period (with 
more than 50 e- foldings) since H^^l^ then corresponds 
to an unobservably large length scale. Therefore, one 
can expect three spectral regimes, {kph/H)i > af/ui 
(always subhubble), af/at > {kph/H)i > [ai/af)''/'^ 
(intermediate) and {kph/H)i < (ai/aj)''/^ (unaffected). 
Translating this to the highest affected mass scale 
involved we estimate Mmax ^ IO^^Mq {a-f /at)^''^'^ . 
Above this scale the spectrum of density perturbations 
is given by the primordial spectrum of density pertur- 
bations, e.g. a nearly scale invariant spectrum. Note 
that we do not make a statement about the detailed 
evolution of perturbations above or below these two 
scales at this point, we only stress that a cosmologically 
interesting mass scale appears for a short period of 
infiation that could lead to observable consequences. For 
a fully consistent treatment of perturbations one needs 
to take into account the dynamics of the chiral phase 
transition in a detailed model and try to estimate the 
effects of reheating on the amplitude of perturbations. 
For cold dark matter the dilution of the energy and 
number densities leads to the possibility of a matter 
dominated phase before the inflationary phase since the 
dark matter energy density after reheating is basically 
fixed by the present day value. Consequently the dark 
matter density before inflation is larger by the same 
factor 9^ as the baryon density. For 6* > 10'^ a matter 
dominated phase is present before the QCD phase transi- 
tion and QCD inflation is naturally limited to a length of 

.1/3 / \4/3/__\l/3 
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The highest affected dark matter mass scale would 
be then Mmax ~ (10^ - 10^)Mq. One can put a 
general upper limit on the amount of entropy that 
is released by demanding that the initial baryon 
asymmetry is at most of order one, implying that 
^max ^ (l/'7B(a_B))^/'^ ~ 1200 (a complete spectrum of 
primordial fluctuations would require 9 > 10^"). 
We note that for non-relativistic decoupling of dark 
matter the weak interaction cross section will no longer 
give the right amount of dark matter today, the dark 
matter annihilation cross section has to be much smaller, 
i.e. CT^m™'' ~ 9-^ a""'"'' as Qdm oc I/ct^J™*'' (we ignore 
logarithmic dependencies on the dark matter mass). 
This gives the interesting prospect that the little infla- 
tion can be probed by the LHC since the discovery of 



a standard weakly interacting massive particle like the 
neutralino would exclude the scenario. 
For thermally decoupled ultra-relativistic par- 
ticles the ordinary temperature relation to the 
radiation background is changed after inflation 

T = TDM^(<7|//(T^i5ec)/ff|//(r)}'^'. General- 
izing the mass limit found in l25|| one arrives at 
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This 



allows for a much higher mass of a thermal relic particle 
without the need for a large number of additional 
effective decoupling degrees of freedom beyond those of 
the standard model. 

A vanishing speed of sound during a first order phase 
transition can also lead to the formation of primordial 
black holes (PBH) for a small fraction of Hubble volumes 
that are sufficiently overdense |26| . The mass spectrum 
of these PBH will be strongly peaked around \Mq which 
corresponds to the total (not just the dark matter) 
energy density inside the Hubble volume at the phase 
transition. The abundance of PBH depends on the 
spectral index and amplitude of the density fiuctuation 
spectrum that can differ significantly around the Hubble 
scale in the presented scenario as discussed above. 
Lumps of quark matter or small quarks stars could be 
also produced but only with M ~ 1Q~'^Mq as we argue 
that nucleation starts after the little inflationary epoch. 
For a first order QCD phase transition with bubble 
nucleation there is a well discussed mechanism for 



producing magnetic fields via bubble collisions [27 1. 
Since the baryon number is carried by massless quarks 
and massive nucleons in the respective phases the baryon 
number will tend to concentrate in the quark phase, at 
least close to the phase boundary [23|- Because of their 
finite masses the muon and the strange quark are already 
slightly suppressed at the critical temperature Tc which 
leads to a charge dipole layer at the phase boundary. 
The resulting net positive charge density is = (ieuB 
with /? - 10^2 _ iQ-i foj. a small t^b and /3 = 0.2 for 
our case. Using the estimates of ref.ji^] we arrive at 
magnetic fields of strength Bqcd = 10* — 10^°G for 
low baryon asymmetry although MHD turbulence may 
readily amplify the initial fields to the equipartition 
value Bf.q = ^SirT'^v'j and refs. therein), where 

Vf is the fluid velocity. In the little inflation scenario 
the initial value of the baryon contrast between the 
two phases can be much higher since nucleons will 
be highly suppressed at T - 170McV/6' - 0.2MeV, 
while for a random walk the baryon diffusion length 
i~diff OC ^1 y/ns + n-B ^ 4//m 9'^^'^ ^ 10 cm is larger 
because ns and are reduced by a factor of 9^. 
Altogether one can expect that the magnetic field B 
will easily reach an equipartition value of B^q ~ lO^^G, 
where Vf ^ 1 since the released latent heat is much 
larger than the thermal energy. 
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The presently observed galactic and extragalactic mag- 
netic fields have strength B^''^ = 0.1 — 1/xG, but the re- 
quired seed fields for an effective galactic dynamo mech- 
anism on scales of 0.1 Mpc are strongly model and pa- 
rameter dependent and vary over many orders of mag- 
nitude 10-30G < ^'^ < 10-i°G (see ^ and ref- 
erences therein). In ref. [30[ it was argued that for a 
causal production mechanism the spectrum of the gen- 
erated magnetic field must be very blue for uncorrelated 
superhorizon scales, i.e. oc A~" with n > 2. There- 
fore B'^'^''' can be strongly limited by the allowed ad- 
ditional energy density at BBN [3l|. We find that the 
produced initial field corresponds to B^'^^^j^^ < 10-22G 
which translates to a mean field at the QCD scale of at 
most Bqcd = 5-lO^^G which is consistent with the above 
estimates. In [slj it was found that an inverse cascade 
mechanism due to a non- vanishing helicity of the primor- 
dial magnetic field (as one can expect in the presented 
scenario) is able to successfully seed large scale magnetic 
fields fields at the QCD phase transition. 

In a first order phase transition nucleation can produce 
gravitational waves due to bubble collisions and hydro- 
dynamic turbulence as found by [ssf . For a nucleation 
rate of F oc exp(i/T) the peak frequency of the spec- 
trum corresponds to a present day frequency of i/^ ~ 



4.0 • 10"**Hz 



1/6 



peak ^ 

where T* is 



i5oMeVy V 50 

the reheating temperature. With the above estimates 
one arrives at a peak strain amplitude ^c(t'^afc) = 4.7 
^2 / isoMeV A ^ ^^/^ 



10-15 I 
lisions. 
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due to bubble col- 



The kinetic energy of the colliding bubbles is 
also partially converted to turbulent bulk motion of the 
plasma stirring gravitational waves at a slightly lower fre- 
quency fjg^jfc — 0.3 i^p^cik with a higher peak amplitude 
^ci^leak) - '^■^^c{Vpeak) ^^"^ ^ strougly first ordcr phase 



transition j34| . The approximate shape of the strain am- 



plitude spectrum is then hc{v) cx j/^/^ for v < H (uncor- 
related white noise) and hc{v) cx i^^'" for v > v^^^f^ where 
the spectral index m is at most 2 but could easily be close 
to 1 or even lower due to multi bubble collisions [35| . Pul- 
sar timing already limits nucleation with the presently 
available data to t/H~^ < 0.12 which will improve to 
r/iJ-i < 0.06 for the full data of the Parkes Pulsar Tim- 
ing Array project Q. The planned Square Kilometer Ar- 
ray will be about four orders of magnitude more sensitive 
in figu, (v) which corresponds to one order of magni- 
tude improvement for the bound on t/H^^. Detection 
via the Laser Interferometer Space Antenna (LISA) could 
also be possible if the high frequency tail of the spectrum 
has a spectral index m < 1.4 and t/H^^ > 10^^. 

We have here only briefly introduced the idea of a little 
inflation at the QCD phase transition and sketched the 
differences from the standard scenario for structure for- 
mation, dark matter properties, magnetic flelds and grav- 
itational wave production. The main assumptions are 



a high initial baryon asymmetry before the QCD phase 
transition and the existence of a quasi stable QCD vac- 
uum condensate that dominates the energy budget for a 
short period. Especially the impact on structure forma- 
tion in this approach seems to be rather interesting but 
requires a more thorough field theoretical approach. 
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